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ABSTRACT 

Wc have observed the ultracompact HII region G29.96— 0.02 in the near 
infrared J, H, and K bands and in the Br7 hne. By comparison with radio 
observations, we determine that the extinction to the nebula is Apc = 2.14 
with a 3(7 uncertainty of 0.25. We identify the ionizing star and determine 
its intrinsic K magnitude. The star does not have an infrared excess and 
so appears to be no longer accreting. The K magnitude and the bolometric 
luminosity allow us to place limits on the location of the ionizing star in the 
HR diagram. The 3a upper limit on the effective temperature of the ionizing 
star is 42 500 K. We favor a luminosity appropriate for star with a mass 
in excess of about 60 Mq. The limit on the temperature and luminosity 
exclude stars on the ZAMS and stars within 10^ yr of the ZAMS. Since the 
age of the UC HII region is estimated to be only about 10^ yr, we suggest 
that this is direct evidence that the stellar birth line for massive stars at 
twice solar metallicity must be significantly redder than the ZAMS. 

1. Introduction 

Our understanding of the formation of massive stars is primitive. In this paper, we 
show how observations of massive stars in ultracompact (UC) HII regions can provide 
information about their birth. UC HII regions are found around young massive stars 
still embedded in their natal molecular clouds. Statistically, the objects classified as 
UC HII regions by Wood & Churchwell (1989a, 1989b) appear to last for about 10-20% 
of the life of an O star or a few x 10^ yr. Massive stellar outflows (Harvey & Forveille 
1988; Garden & Carlstrom 1992; Hunter ct al. 1995; Shepherd & Churchwell 1996a, 
1996b) seem to represent a still earlier and more poorly understood phase, possibly 
during which the massive star is still accreting. 

Previously, the ionizing stars in UC HII regions have been studied through their indirect 
effects on the surrounding gas and dust. These studies include determinations of 
the bolometric luminosities (Chini et al. 1986; Wood & Churchwell 1989a), effective 
ionizing fluxes (Wood & Churchwell 1989a; Kurtz, Churchwell, Sz Wood 1994), and 
effective spectral hardnesses (Doherty et al. 1994). All of these methods have significant 
drawbacks. The bolometric luminosities are upper limits because they arc determined 
largely from IRAS observations and most likely include contributions from luminous 
sources other than the ionizing star of the UC HII region. The effective ionizing flux 
and the effective spectral hardness are very difficult to relate to the star because of 
the effects of diist within the UC HII region and large uncertainties in stellar ionizing 
continua (see the discussion in Najarro et al. 1996). 
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Instead, we have chosen to observe the UC HII region G29.96— 0.02 in the near infrared 
J, -H, and K bands and in the 2.17;um Br7 hne. These observations ahow direct 
measurements of the properties of the ionizing star. G29.96— 0.02 has been previously 
imaged in the near infrared by Van Buren (1993), Megeath (1993), and Fey et al. (1995). 
These authors extracted a great deal of qualitative information from their observations 
and were able to make a number of suggestions. Our work builds on their work; we are 
able to confirm a number of their suggestions and derive important new information 
on the ionizing star and the nebula. 

Fey et al. (1995) considered a champagne flow model for G29.96— 0.02 and estimated 
an age for the UC HII region of about 10^ yr from the extent of the extended radio 
continuum emission. Mac Low et al. (1991) and Van Buren & Mac Low (1992) modeled 
029.96— 0.02 as a bow shock around a star moving at 5-20 kms~^ at an angle of about 
135° to our line of sight. However, the O star is closely associated with a young cluster 
(Fey et al. 1995 and below) and a hot core apparently containing a protostar (Cesaroni 
et al. 1994); this suggests that it could have moved at most 10 arcsec and so, again, has 
an age of less than 10^ yr. Lumsden & Hoare (1996) found that a bow shock is a poor 
fit to their Br7 velocity field and prefer instead a model similar to the one suggested 
by Fey et al. (1995). Mass- loaded models (Hollenbach et al. 1994; Dyson 1994; Lizano 
& Canto 1995; Dyson, Williams, &; Redman 1995; Redman, Williams, & Dyson 1996; 
Williams, Dyson, & Redman 1996; Lizano et al. 1996) can delay the expansion of an 
UC HII region by about 10^ yr. If the central star was in an outflow/accretion phase 
before the UC HII region was created, this phase is likely the have lasted less than 
10^ yr, based on comparison with other O stars currently undergoing or just finishing 
molecular outflow activity (Cabrit & Bertout 1992; Shepherd & Churchwell 1996b; 
Acord, Walmsley, &; Churchwell 1997). All of these estimates and considerations point 
to an age of order 10^ yr for the UC HII region. 

We adopt a heliocentric velocity of 95 kms"-*^, the average global radio recombination 
line velocity (Afflerbach et al. 1994). Since G29.96— 0.02 lies close to a tangent point, 
its distance is uncertain, and, unfortunately, the formaldehyde absorption line data 
of Downes et al. (1990) do not resolve the distance ambiguity. The Galactic rotation 
curve and error bars given by Burton (1988) yield a distance of between 5 kpc and 10 
kpc. 

Afflerbach, Churchwell, & Werner (1996) have determined that the abundances of O, 
N, and S in G29.96— 0.02 are about twice those observed in the local ISM. This is 
consistent with its location towards the center of the galaxy. 

The organization of this paper is as follows. In section 2 we present our new infrared 
observations. In section 3 we summarize the existing radio observations that have been 
kindly made available to us by a number of our colleagues. In section 4 we present the 
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astrometric calibration of the location of the ionizing star. In section 5 we discuss the 
bolometric luminosity of G29.96— 0.02. In section 6 we investigate the nebula using 
the Br7 and radio images and derive the extinction to the nebula. In section 7 we 
investigate the properties of the ionizing star and its associated stellar cluster using 
the broad band J, H, and K images. In section 8 we summarize the principal results 
of this work and discuss some of their implications. 

2. Infrared Data 
2.1. Observations 

We observed G29.96-0.02 on the night of 1995 May 21 and June 24 UT with the 
OSIRIS near infrared camera (Depoy et al. 1993) on the Perkins 1.8 meter telescope at 
Lowell Observatory. Light cirrus was present during the first night. The second night 
was photometric. The seeing was 1.8 arcsec FWHM on the first night and 2.0 arcsec 
FWHM on the second night. 

OSIRIS has a 256x256 NICMOS-3 detector. The re-imaging optics gave a pixel scale 
of 0.62 arcsec and a field of view of 150 arcsec. Images were taken through J, i7, and 
K broad band filters and 2.12, 2.14, and 2.17 /im narrow band filters. The 2.12 ^um 
filter includes emission from the H2 = 0-1 S(l) line and the 2.17 ;um filter includes 
emission from the Br7 line; the 2.14 /xm filter excludes strong emission lines. 

We took several exposures through each filter, dithering the target onto each of the 
detector quadrants. Total exposure times in each of the broad band filters were 120 
seconds on the first night and 240 seconds on the second night. Total exposure times 
in each of the narrow band filters were 480 seconds on the first night and 240 seconds 
on the second night. 

2.2. Reduction 

The data were reduced using custom programs within the VISTA image reduction 
system. Each image was corrected for non-linear response by subtracting a bias image 
and dividing each pixel by its estimated relative response. We have found that the 
non-linearity can be adequately modeled as a quadratic function and the principal 
variation in non-linearity in OSIRIS is between quadrants, so we determined and 
applied different quadratic coefficients for each quadrant. Each image was then dark 
subtracted, fiattened, and masked to eliminate bad pixels. A single sky value was 
estimated for each image by fitting a quadratic to the peak of the intensity histogram 
and solving for the turning point. That value was then subtracted from the image. 
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The images were brought into a common registration by measurements of bright stars, 
shifted by an integral number of pixels, and co-added. This process did not noticeably 
degrade the seeing. 

We corrected the 2.12 and 2.17 fim images for continuum emission by subtracting the 
2.14 yum image after scaling by the ratio of the mean brightnesses of several stars away 
from the UC HII region. We did not detect any line emission in the 2.12 fim image 
and will not consider it further. 

2.3. Broad Band Calibration 

The broad band data from the second nights data were calibrated by observations of 
the two UKIRT faint standards 25 and 28 (Casah & Hawarden 1992). The UKIRT faint 
standards contain very little color information, so we did not attempt to solve a color 
term and have left the magnitudes on the natural system. This is less of a drawback 
than it seems, as our sources are so heavily reddened that color terms derived from 
observations of unreddened stars would be of limited use. The RMS residuals in the 
fits for the zero points and extinction terms were about 3%. 

2.4. Narrow Band Calibration 

The narrow band data from the second night was calibrated by observations of the 
sdO optical spectrophotometric standard BD +28 4211. This star is not an infrared 
spectrophotometric standard but we expect that its spectrum will be featureless in 
the K window. BD +28 4211 has a fainter red companion 3 arcsec to the NNW 
(Massey & Gronwall 1990; Thejll, Theissen, & Jimenez 1994). We measured a difference 
between the fainter and brighter stars of AK ~ 1.5. Before performing photometry, 
we removed the companion by shifting, scaling, and subtracting a suitably centered 
star. We measured a K magnitude of 11.56 ± 0.04 for BD +28 4211. The absolute 
calibration of Bessell & Brett (1988) gives a flux of 9.62 x 10^^^ ergs"^ cm~^ at 
2.19 //m with a la uncertainty of about 4%. Assuming a A^^ spectrum, this gives a 
flux of 1.01 X 10~^^ ergs~^ cm~^ at 2.164 ^um, the peak of the 2.17 micron filter. 

2.5. Infrared Images 

Figure || presents images in J, H, K, and Br7. The arc-shaped nebula at the center 
of each image is the UC HII region G29.96— 0.02. The bright point source within the 
nebula in the J, H, and K images is the ionizing star candidate previously seen by 
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Van Buren (1993), Megeath (1993), and Fey et al. (1995). The position of this source 
is marked with a white cross in the Br7 image. 

3. Radio Data 

G29.96— 0.02 has been extensively studied in the radio continuum, radio recombination 
hues, and molecular lines. We are grateful to a number of our colleagues for providing 
the radio images listed in Table |^ in electronic form. Table |l] lists for each image the 
authors, frequency, telescope, beam FWHM, and largest angular size imaged. The VLA 
images are all of radio continuum emission. The OVRO image is the total intensity in 
the C^^O J = 1 — line. We assume that the 3a errors in the absolute flux calibration 
of the radio continuum images are about 20%. The errors in the C^^O image are larger 
but do not affect our analysis. 

4. Astrometry 

We determined the position of the ionizing star candidate (#3 below) by cross- 
correlating versions of the Br7 image with the astrometrically-calibrated 1.3 and 2 
cm radio continuum images. Before cross-correlating, we smoothed each pair of images 
to the same resolution. The ionizing star candidate is located at aigso = 18 43 27.214 
and (^1950 = —02 42 35.9 in the radio reference frame. The determinations from the 
1.3 cm and 2 cm images differed by only 0.2 arcsec. The distance between the ionizing 
star and the peak of the 2 cm emission is 2.0 arcsec. At 5.0 kpc this corresponds to 
10000 AU. 

5. The Bolometric Luminosity 

Table |2| lists measurements of the continuum spectral energy of G29.96— 0.02. The 
values at 21 cm, 6 cm, and 1.3 cm were measured by us from the VLA maps obtained 
by Claussen & Hofner (1996), Afflerbach et al. (1994), and Cesaroni et al. (1994). The 
three measurements at 1.3 mm were obtained with different effective beam sizes. The 
measurements at 790 //m and 350 fim from Hunter (1997) may be missing extended 
emission. To obtain total magnitudes of J= 12.5, H = 10.1, and K = 7.9 for the nebula 
we measured the total magnitudes within a 25 arcsec box centered on the ionizing star 
after subtracting scaled point spread functions of all stars except the ionizing star. (We 
describe our photometry below.) We converted these to Jy using the flux calibration 
of Besseh & Brett (1988). Our measurements of AJ 1.9, AH 1.8, and AK w 2.4 
between the flux for the ionizing star and the total flux show that nebular emission and 
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scattering dominate over direct light from the star even at these short wavelengths. 

Figure ^ shows the spectral energy distribution plotted as log vF,^ against log u. The 
spectral energy distribution is characteristic of UC HII regions: optically thick free-free 
emission longward of a few cm, optically thin free-free emission between a few cm to 
about a few mm, and very strong thermal dust emission from about 1 mm to about 
10 nm. 

We can estimate the bolometric luminosity by integrating over the spectral energy 
distribution. As expected, extrapolation of the cm wavelength free-free emission to 
1300 /im (the dashed line in Figure |2|) indicates that the 1300 fim emission arises from 
hot dust. We fitted a modified Planck law F,y oc v'^By{T) to the 1300 ^m and 100 //m 
points, finding T ~ 26 K, and interpolated elsewhere. This adopted spectrum is shown 
by a solid line in Figure |2|. We also considered the optically thick model fit by Hunter 
(1997); this gives a flux within 5% of ours. We derive an apparent bolometric flux of 
7.00 X 10"'' ergs~^cm~^ which corresponds to an apparent bolometric magnitude of 
"^boi = 3.90 (where Mboi,© = 4.75). Since this depends heavily on the IRAS 100 jiva 
measurement, we assume a 3cr uncertainty of 30%. 

Since this bolometric luminosity is integrated over an arcmin sized region, it is likely to 
include contributions from other sources and so is a priori only an upper limit for the 
bolometric luminosity of the ionizing star. However, there are no other bright sources 
in either the 1300 /xm map of Mooney et al. (1995) or the mid infrared images of Ball 
et al. (1996). We do know that the IRAS beam includes both the hot core to the west 
of the UC HII region and a cluster of B stars apparently associated with the ionizing 
star (see below). Cesaroni et al. (1994) estimate that the source powering the core 
contributes only about 1/10 of the flux. Hofner et al. (1997) estimate that at 2.7 mm 
only 150 it 50 mJy of the total 2.73 Jy arises from the hot core. Much of the total 
emission will be free-free emission from the UC HII region, but since the hot core does 
not appear in cm wavelength continuum maps we attribute its emission solely to dust. 
We estimate the total dust emission of the region at 2.7 mm to be about 4 Jy by 
extrapolating the 1.3 mm measurement of 15.6 Jy (Mooney et al. 1995) assuming a v"^ 
spectrum. Thus, the hot core contributes less than 10% of the 2.7 mm dust emission. 
We expect that the hot core has a spectrum that is no hotter than the UC HII region and 
so this fraction is an upper limit to its total contribution to the bolometric luminosity. 
Furthermore, as we shall see below, the cluster appears to be very much less luminous 
than the ionizing star. This suggests that the true bolometric luminosity of the ionizing 
star is likely to be quite close to the total bolometric luminosity of the region. We will 
conservatively assume that at least half of the measured luminosity must come from 
the ionizing star, that is 3.90 < m-boi ^ 4.65 with a 3(T uncertainty on the limits of 0.3. 
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6. The Nebula 

6.1. Emission Measure and Electron Temperature 

In order to predict the intrinsic flux in Br7, we need to know both the emission measure 
and the electron temperature in the nebula. We can use the 2 cm, 6 cm, and 21 cm 
radio continuum images to provide this information. The familiar expression for the 
brightness temperature Tb at a given frequency v is 

Tb = re(l-e-"). (1) 

The optical depth is given by 

r = 8.235 X \^-'^aT^"^-^^(y jQ'RzY'^-^iE jciTi'^ pc) (2) 

where E is the emission measure and a is a weak function of v and (Mezger & 
Henderson 1967). The peak synthesized beam brightness temperature measured in the 
2 cm image is about 1500 K. Anticipating our result that Tg is about 6000 K and 
assuming that there is no unresolved structure in the 2 cm image, equation 1 implies 
that the peak optical depths at 2 cm, 6 cm, and 21 cm are about 0.3, 3, and 40. Thus, 
the 6 cm and 21 cm images provide information on different regions of the nebula. 
Since the dependence of o on Te is so weak, we again anticipate our result and use its 
value at 6000 K. 

We smoothed the 2 cm image to match the resolution of the 6 cm and 21 cm images and 
then calculated Tg and E iteratively. We started with an initial guess for Tg of 6000 K. 
We then used the 2 cm image and equations 1 and 2 to predict the 6 cm or 21 cm image. 
Half of the ratio between the predicted and actual 6 cm or 21 cm image was used at 
each point in the image as a correction factor for the value of Tg. We repeated the 
prediction/correction step until Tg no longer changed significantly. Images of Tb and 
Tg are shown in Figure |^. The effects of increasing optical depth at longer wavelength 
are apparent as increases in Tb and a severe reduction in contrast in the 21 cm Tb 
image compared to the 2 cm Tb image. In the optically thick part of the nebula, Tg is 
determined largely by the brightness temperature in the longer wavelength image and 
so has a systematic 3(T error of about 20%. In the optically thin parts of the nebula, 
it is difficult to separate the product of Tg and E and the errors in Tg are much larger. 
Random errors from imperfect imaging are more difficult to judge. The hot regions 
along the NW and SW of the arc in Figure |3|c and the 'bumps' in Figure ^f are almost 
certainly artifacts. Ignoring these, Tg is between 5000 K and 7000 K over the bright 
part of the nebula, including the arc, and might possibly be about 1000 K cooler in 
the tail to the NE. In our following analysis we will adopt a uniform temperature of 
6000 K and consider temperature variations in our discussion. 
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Afflerbach et al. (1994) examined the nebular properties by a non-LTE analysis of radio 
recombination lines. They derived Tg in the range 6200-8600 K, although with large 
errors. Thus, our measurement of T^. of about 6000 K is in rough agreement. 

6.2. The Extinction To The Nebula 

Figure ^ shows our Br7 image of G29.96— 0.02. The total Br7 flux in Figure Q 
is 8.42 X 10^^^ ergs~^cm~^ with a Scr uncertainty of about 12%. We measure 
2.23 X 10^^^ ergs~^cm~^ in a 5 arcsec aperture centered on the peak of the Br7 
emission, consistent with the flux of 2.49 x lO"^'^ ergs~^cm~^ measured by Doherty 
et al. (1994), and 6.6 x 10~^^ ergs^^ cm^^ in a 17 arcsec aperture, consistent with the 
flux of 7.6 lb 1.0 x 10^^^ ergs~^ cm^^ measured by Herter et al. (1981). 

Figure ^b shows the 2 cm image of Fey et al. (1995) at its full resolution and Figure §c 
shows it smoothed to the 1.8 arcsec FWHM resolution of our Br7 image. As we 
described above, the alignment of the 2 cm image was determined by cross-correlation 
with the Br7 image, after smoothing the 2 cm image to match the resolution of the 
Br7 image. The total 2 cm flux density in Figure ^ is 4.6 Jy. Fey et al. (1995) quote a 
flux of 3.9 Jy for the whole image, but this includes large regions with small, negative 
flux. The similarities between the Br7 and 2 cm images are striking; both show the 
strong arc of emission that lead Wood & Churchwell (1989a) to classify G29.96— 0.02 as 
a cometary UC HII region and both show emission behind the arc and in extensions to 
the south east and the north. The importance of this is that it demonstrates that the 
emission occurs at similar optical depths and so the extended emission noted by Fey 
et al. (1995) is directly associated with the UC HII region rather than being a chance 
superposition of foreground or background emission. 

We now derive a relation between the 2 cm brightness temperature and the intrinsic 
flux in Br7, with the purpose of comparing the intrinsic Br7 flux to the observed Br7 
flux to derive the apparent extinction at Br7. Unless otherwise noted, all units are cgs. 
From equations 1 and 2, the emission measure E in cm~^ is 

E = ^.na~,^^T^-^'vl-^ In (^r^) • (3) 

The Br7 flux 5'Br7 from a region of solid angle 17 is 

5Br7 = G.9huBr^af,^^E. (4) 

The factor of 0.9 arises because we have assumed that the nebula consists of 10% singly 
ionized helium by number; the ionized helium will contribute electrons and ions to the 
radio continuum emission but only electrons to the Br7 emission. Over the range of 
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temperatures of interest here, Hummer &: Storey (1987) give 

af,^ = 6.48 X 10-"re-i-°^ (5) 

We ignore variations of a^^^ with density as they amount to less than 1% over the 
range of interest here. Combining equations 3, 4, and 5 and replacing a2cm and i'2cm 
gives 

To understand the errors it is worth making the approximation that Tb <C T^. The 
peak synthesized beam brightness temperature in the 2 cm image is only 1500 K, so 
this is a good approximation to first order. We then have 

^Br^ ~ o.6i2^re-o-^^rb. (7) 

It can be seen that the 20% Scr systematic uncertainty in Tb and Tg and the 12% 3a 
uncertainty in the observed flux in Br7 translate to a 27% 3a systematic uncertainty 
in the extinction. 

We used equation 6 to estimate the intrinsic flux of Br7 at each point from the 2 cm 
image. The total intrinsic flux in Figure § is expected to be 6.72 x 10~^^ ergs~^ cm~^ 
and so the apparent extinction at Br7 is 2.26 mag. This extinction is an average over 
the whole nebula and we quote it only for illustrative purposes; we compute images of 
apparent extinction below. 

The apparent extinction to the nebula is not uniform. Figure ^d shows an image of the 
apparent extinction. The apparent extinction is higher in the arc and towards the end 
of the extension to the south east, reaching about 2.6 mag. In the region behind the arc, 
in the vicinity of the ionizing star, the apparent extinction drops to about 2.20 mag, 
although there is some noise from continuum subtraction. The peak-to-valley contrast 
in the main part of the nebula is about 0.6 mag or a factor of about 1.7. It seems 
significant that the two regions of highest apparent extinction correspond to the two 
brightest regions. We can think of six possible explanations for the correlation of higher 
apparent extinction with surface brightness: an error in the 2 cm map, the presense 
of dense ionized regions in the nebula, variations in electron temperature within the 
nebula, variations in external extinction, a molecular sheath, and internal extinction. 
We consider these in turn. 

(i) Errors in the 2 cm image. The kind of error in the 2 cm image that would explain 
the structure would have the characteristic of placing almost twice as much flux in 
small scale structures (the arc) than large scale structures (the tail). Since the 2 cm 
image is constructed from images in three VLA configurations, such an error might have 
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occurred if there were relative flux errors between the observations. To investigate this, 
we compared the 2 cm image to the 1.3 cm image. The peak optical depth in the 2 cm 
image is 0.3 and so optical depth effects will be important at the 10% level in the arc. 
Correcting both images for optical depth effects is complicated, since they have different 
intrinsic resolutions. Therefore, we compare the real 1.3 cm image to the prediction 
from the 2 cm image, obtained by assuming Tg = 6000 K, using equations 1 and 2 
to account for optical depth and frequency effects, and smoothing from the resolution 
of the 2 cm image to that of the 1.3 cm image. Figure ^a and b show the observed 
and predicted 1.3 cm images. The predicted image has only 87% of the flux in the 
observed image, but this is not unexpected. Although there are significant variations 
in the ratio, they are only at the 10% level and are not in the sense required to explain 
the additional apparent extinction (the arc does not appear brighter in the 2 cm image 
than the 1.3 cm image). It would appear, then, that errors in the 2 cm image cannot 
explain the variations in apparent extinction. 

(ii) Dense Ionized Regions. Our procedure above accounts for optical depth effects as 
long as there is no unresolved structure. It is possible, though, that unresolved dense 
ionized regions exist in the nebula. These would contribute strongly to Br7 but might 
well be optically thick at 2 cm and contribute only weakly. However, in order to mimic 
the variations in the apparent extinction image, these clumps would have to avoid the 
brightest regions in the nebula. This seems somewhat unlikely. 

(iii) Variations in Electron Temperature. As equation 7 shows, the Br7 flux depends 
on the electron temperature roughly as T^""'^^. Thus, we could explain the factor 
of 1.7 peak-to- valley contrast in the ratio by a factor of 2 peak-to- valley contrast in 
the temperature, in the sense that the brighter regions are cooler. However, in our 
investigation of Te above we discovered that this was not the case: the region around 
the arc is certainly not cooler than the tail and is perhaps 1000 K or 15% hotter. 

(iv) Variations in External Extinction. At first glance, variable external extinction to 
the nebula seems a likely candidate. Figure ^ shows that G29.96— 0.02 is superposed 
on a molecular clump with a dense core in front of the arc. However, the contours 
of molecular gas do not match the details of the variation in extinction: they would 
suggest a smooth gradient across the nebula rather than the strong correlation with 
nebular emission that we observe. 

(v) A Molecular Sheath. The bow shock model predicts a compressed molecular sheath 
around the arc. Van Buren and Mac Low (1992) used the radio recombination line 
observation of Wood &: Churchwell (1991) to show that the velocity field implied a 
viewing angle of about 135° in the bow shock model, i.e., the head is more distance 
than the tail. Nevertheless, it is possible that we are viewing the edge-brightened arc 
at an oblique angle through such a sheath. 
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We can estimate the column density through such a sheath by adopting Ay/Nu = 
1.0 X 10^^^ magcm"^ and A-Q^^/Ay = 0.125 (for the Ry = 5.0 extinction law of 
Cardelli, Clayton, & Mathis 1989). Our value of Ay/N]i is twice the local value to 
account for the higher metallicity in G29.96— 0.02. We take the local average from 
Bohlin, Savage, & Drake (1978), although Cardelh, Clayton, & Mathis (1989) give 
reasons for believing that the value may be lower in molecular clouds. The excess A-qt-^ 
of about 0.4 mag corresponds to Ny{^ of about 1.6 x 10^^ cm~^. This column is similar 
to the predictions of A^H2 ~ 5 x 10^^ cm~^ for the column density through the shell 
in the vicinity of the stagnation point (Van Buren et al. 1990; Van Buren & Mac Low 
1992), although the comparison is simplistic because we ignore the effects of geometry. 

(vi) Internal Extinction. Internal extinction would provide a natural explanation of the 
correlation between excess apparent extinction and nebular emission. We can eliminate 
radiation transfer effects in hydrogen by noting that the peak emission measure in the 
2 cm image is 1.4 x 10*^ cm~^ pc and thus the line center optical depth in Br7 is less 
than 3 x 10~^ (Hummer & Storey 1987). However, dust may provide significant opacity. 
We can investigate this by assuming that dust is evenly mixed with a uniform density 
ionized gas in the nebula. We ignore the effects of scattering, which may be a better 
approximation than might be supposed, since the arc is most likely to be a sheet in 
which scattering will effectively remove photons from the line of sight. In this model, 
the internal extinction is 

J^^,^ = 0.5 (^Br^Mv) {Ay/Nn) Nn (8) 

The factor of 0.5 arises because the effective optical depth of an evenly mixed slab of 
emitting matter is — In Jq e~'^^ dl = — ln((l — e~'^)/r) = r/2 + O(r^). If the nebula 
is uniform, then A^h = EM/n^ and we can derive A^h from the 2 cm image and the 
electron density of 6 x 10^ cm~^ (the average of the values from regions A, B, and W from 
AfRerbach et al. 1994). The peak emission measure in the full-resolution 2 cm image is 
1.4 x 10^ cm~^ pc and so this model predicts a peak effective internal extinction of 0.41 
magnitudes. However, the extinction must be weighted by the emission and smoothed 
before it can be directly compared to the Br7 map. When this is done, the peak effective 
internal extinction drops to 0.27. The emission-weighted value of A!^^^ is shown in 
Figure |^. Figure |^b shows A-^^^ — Agj,^, the apparent extinction after subtracting 
the model internal extinction. The pattern of higher apparent extinction is markedly 
reduced in the vicinity of the arc but is not eliminated. Figure 0c shows A-q^^ — 2A'^^^, 
the apparent extinction after subtracting twice the model internal extinction. This 
almost completely removes the traces of non-uniform extinction over the majority of 
the nebula. A number of things could act to raise the internal extinction above our 
model. We could have underestimated Ay/N^i by a factor of two, but this seems 
unlikely (Cardelli, Clayton, & Mathis 1989). It seems more likely that the nebula is is 
not uniform in density but includes relatively dense ionized regions with Ue ~ lO^cm"'^ 
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that contribute to both the emission and extinction along with either lower density 
ionized regions or neutral regions that only contribute to the extinction. 

Ultimately, our lack of an alternative explanation and the striking correlation between 
the variations in extinction and the brightness of the nebula leads us to the conclusion 
that we must be seeing the effects of either a molecular sheath around a bow shock 
or of internal extinction. That our model for internal extinction fails to predict the 
magnitude of the extinction by a factor of two suggests that there are lower density or 
neutral inclusions in the arc of the nebula with about as much column density as the 
dense ionized gas. This might naturally occur in models in which UC HII regions ablate 
dense inclusions (Dyson 1994; Lizano & Canto 1995; Dyson, Williams, & Redman 1995; 
Redman, Williams, & Dyson 1996; Williams, Dyson, & Redman 1996; Lizano ct al. 
1996) or in which the arc arises from a photoevaporation front being driven into an 
inhomogeneous medium. 

The location of the star slightly off the arc means that, regardless of the origin of the 
additional apparent extinction associated with the arc, the extinction to the ionizing 
star at Br7 is 2.20 mag with a 3a uncertainty of 0.25 mag. 

7. Stellar Properties 

7.1. Near Infrared Photometry 

Stellar photometry in the region of the UC HII region is hindered by the strong and 
variable nebular background, the degree of crowding, and the extreme colors of some 
of the stars which cause them to be undetected in one or more of the three bands. We 
adopted methods to minimize these problems. 

First, we performed stellar photometry on the J, H, and -flT images using Jon Holtzman's 
very heavily modified version of DAOPHOT (Stetson 1987). One of the modifications 
allowed us to fit stars in the three bands simultaneously, fixing the relative positions 

of stars within groups but allowing the offsets between frames to vary. We carefully 
added stars to each group by hand and did not allow the software to delete any. Wc 
restricted ourselves to the region of common overlap between our dithered exposures, 
to avoid problems associated with PSF variations (McCaughrcan 1993). 

Second, wc scaled and subtracted the Br7 line image from each of the broad band 
images before performing photometry. The scaling was determined subjectively. This 
procedure works well, but it is not perfect. Clearly, the ratio of broad band nebular 
emission to Br7 emission is not constant, presumably because of the presense of a 
component such as scattered light or dust emission that is not proportional to the 
emission measure and because of variable extinction. The effect of this procedure can 
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be seen in Figure y which shows images of the nebula in J, H, and K both before and 
after suppression of the nebular emission. 

Our photometry depended to some degree on the scaling factors applied to the Br7 
line image. We determined this additional uncertainty by performing photometry on 
images for which the scaling factors had been changed by ±10%. Since these changes 
produced images that were quite noticeably under- or over- subtracted, we assumed 
that they correspond to Scr errors. For the ionizing star these errors amounted to a 
0.01-0.02 mag contribution to the lex error but for the faintest stars they amounted to 
as much as a 0.3 mag contribution to the lex error. 

Our first night's data were taken under better seeing conditions than the second's, but 
the second night's data were photometric. Accordingly, we derived relative magnitudes 
from the first night's data and fixed the zero point using aperture photometry of 
relatively isolated stars in the second night's data. 

Table |3| gives astrometry and photometry for all stars in a 1 arcmin box centered on 
the ionizing star that have photometric errors of less than 0.5 in K and 1.0 in both the 
J—H and H—K colors. The errors shown are la. Table ^ is ordered by increasing K 
magnitude and includes the offsets (Ax and Ay in arcsec) from the ionizing star. The 
ionizing star is ^^3 and its close neighbor to the WSW is ^^4. 

We measure K= 10.36 ±0.04, J-H= 2.40 ±0.05, and H-K= 1.59 ±0.05 the ionizing 
star. Lumsden & Hoare (1996) report Kn = 11.2, J—H = 2.4, and H — Kn = 1.6 
with "typical errors of 0.3-0.5 mag". Thus, their colors are in good agreement with 
ours but their magnitudes are about 0.8 mag fainter. We assessed the reliability of our 
broad band photometry in the following manner. Our Br7 flux is in good agreement 
with other measurements, suggesting that our narrow band calibration is essentially 
correct. We measured 2.17 /xm fluxes for the isolated stars #1, #2, and #6 using simple 
aperture photometry. These fluxes agreed with our K magnitudes to within 0.02 mag. 
This suggests that the our K photometry of isolated stars is essentially correct. We 
do not believe we could have made a 0.8 mag error in the relative magnitudes of star 
#3 to the others; for star #3 to be as faint bs K = 11.2 would require it to be almost 
as faint as star (located 14 arcsec WSW of the ionizing star in Figure ^) but it is 
clearly much brighter. Although we cannot offer an explanation for the disagreement 
between our magnitudes and those of Lumsden & Hoare, these considerations give us 
confidence that ours are more reliable. 

Since the stars were selected by hand, completeness is not well defined. Furthermore, 
the completeness varies considerably between those regions far from the nebula and 
bright stars, where stars as faint as J ~ 17, ~ 16, and ~ 15 have lex errors of 
0.1 and are easy to detect, to those regions close to the nebula and bright stars where 
similarly bright stars might easily have been missed. 
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Figure ^ shows color-color diagrams of stars from Table 1. Only stars with la errors of 
less that 0.15 in K and 0.3 in both colors are shown. Figure ^ shows stars within a 30 
arcsec square centered on the ionizing star and Figure ^3 shows stars outside this square. 
In Figure ^a, the ionizing star is shown with a solid square. The error bars show la 
errors. The dotted lines are the loci of unreddened main sequence and giant stars from 
Bessell &: Brett (1988) and Koornneef (1983). The dashed lines are reddening vectors 
corresponding to Ax oc A^^'^ (lower), A"^ *^ (middle), and A"^''^ (upper) extinction laws 
(Martin & Whittet 1990). We have assumed mean wavelengths of 1.25 ^m, 1.65 ^m, 
and 2.20 ^m for the J, H, and K filters. The length of each vector corresponds to 
^Br7 = 2.20. The quadrilateral delimits the expected colors for a hot star in the 
UC HIT region with Abv-j = 2.20 ± 0.25 or Ak = 2.14 ± 0.25. (The scaling between 
A-Bvy and Ak does not depend significantly on the particular extinction law.) It can 
be seen that most of the stars have colors consistent with moderately reddened normal 
stars. 



7.2. The Ionizing Star 

Van Buren (1993), Megeath (1993), Fey et al. (1995), and Lumsden & Hoare (1996) 
have suggested that star ^^3 is the principal ionizing source of the UC HII region. Star 
#3 has K= 10.36 ±0.04, J-H= 2.40 ±0.05, and H-K= 1.59 ±0.05. These colors are 
in excellent agreement with those predicted for a hot star reddened by the measured 
extinction of the nebula (shown as a quadrilateral in Figure ^a). Of the stars that 
are candidate members of the young cluster (discussed below), #3 is the brightest, 
has no infrared excess, and is situated in the center of the brightest arc of nebular 
emission. This strongly supports the previous suggestions that it is the primary source 
of ionization. There is no evidence for an infrared excess at K. This suggests that any 
remaining disk is optically thin at temperatures of about 1300 K and so the star is no 
longer accreting. 

We can estimate the intrinsic apparent K magnitude rriK of the ionizing star by using 
the extinction of Ak = 2.14 with a 3(T uncertainty of 0.25 and our measurement of 
K = 10.36 with a 2>a uncertainty of 0.12. We find the intrinsic rriK is 8.22 with a 3(7 
uncertainty of 0.28. Since the star does not have an infrared excess, we can estimate the 
intrinsic apparent V magnitude m y from rriK by assuming an intrinsic V— K close to 
—0.90 (Koornneef 1983). We find the intrinsic my is 7.32 with a 2>a uncertainty of 0.28. 
In Section ^ we determined 3.90 < mboi ^ 4.65 with a 3(T uncertainty on the limits of 
0.3. To proceed we need to relate my and mtoi- We adopt the bolometric correction 
BC = JTiboi — rny given by Vacca, Garmany, & Shull (1996). (We use their equations 3 
and 5 and hold the luminosity class constant at V, as the luminosity class makes very 
little difference.) This comes from a fit to a number of non-LTE atmospheres. 
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With our limits on the distance, my, mhol, and my — m-bol we can construct hmits on 
the Lboi and Teg of the star. The regions allowed by la, 2a and 3a errors are shown 
in Figure ^ by thick solid lines. The thin dotted lines show the loci where my = 7.32 
at different distances. Also shown in Figure [l^ are evolutionary tracks and isochrones 
of the Z = 2Zq models of Meynet et al. (1994). The thin solid lines are tracks for 120, 
85, 60, 40, 25, and 2OM0 models, the thick dashed line is the zero-age main-sequence 
(ZAMS), and thin dashed lines are the 1, 2, 3, 4 x 10^ yr isochrones. We only show 
the tracks on the first excursion to the red. The dashed-dotted lines are approximate 
example stellar birth lines. Stellar birth lines have not yet been constructed at Z = 2Zq. 
Instead, we have applied the differences in the Lboi and Teg between the Meynet et al. 
(1994) Z = Zq and Z = 2Zq ZAMS to the Z = Zq 'case 1' (right) and 'case 3' (left) 
models of Bernasconi &: Maeder (1996). (The 'case 1' and 'case 3' models differ in their 
accretion rate.) This difference accounts for the change in atmospheric structure, to 
first order, but fails to account for the increase of the rate of evolution with metallicity 
(Bernasconi 1996, private communication). These birth lines are, therefore, almost 
certainly somewhat too blue. 



It can be seen from Figure |10| that the 1, 2, and 3a upper limits on the effective 
temperature of the ionzing star are 37500, 40 000, and 42 500 K. If the ionizing star 
is a single star or an unequal-mass binary, our observations place a firm lower limit on 
its mass of about 30Mq. If almost all of the bolometric luminosity in the region arises 
from the ionizing star, as we suspect, then this limit is about 6OM0. 

The ZAMS lies entirely outside the region allowed by our observations. Although equal- 
mass binaries are common among O stars (Garmany, Conti, & Massey 1980), a pair of 
binary ZAMS stars are also inconsistent with our observations. For consistency with 
the ZAMS, we require three ZAMS 40Mq stars at the minimum allowed distance and a 
conspiracy of 3a errors. We examined the residual image after subtracting a scaled PSF 
and found no evidence for elongation. We conservatively place an upper limit on the 
separation of the components of such a system as half the FWHM or about 4000 AU 
at 5 kpc. This is much smaller than the size typical for early-type Trapezium systems 
(Abt 1986). The dynamical time for such a compact system is only about 10^ yr, so 
it would dissolve over the estimated age of 10^ yr unless it were formed in a stable, 
hierarchical state. We consider this explanation forced and extremely unlikely. 

The tracks of evolved stars with ages in excess of 10® yr coincide with the region allowed 
by the data, but we reject this possibility because the age of the UC HII region is only 
about 10^ yr. Since the ionizing star appears to be too cool at a given luminosity to be 
a ZAMS star, one might think it could be a pre- main-sequence star evolving onto the 
ZAMS from the Hayashi track at roughly constant Lboi but increasing T^fi- However, 
the modern theory of the formation of massive stars is that pre-main-sequence stars 
evolve almost along the ZAMS with both Lboi and T^q increasing, although they diverge 
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from the ZAMS at the highest masses because of evolutionary effects (Beech &: Mitalas 
1994; Bernasconi &; Maeder 1996). When they finish accreting, the most massive stars 
are somewhat cooler than the ZAMS. This locus is known as the stellar birth line. Our 
observations suggest that the stellar birth line at Z = 2Zq must be sufficiently red that 
the star could evolve from it into the region allowed by the data in only about 10^ yr. 
Such a birth line must be cooler than about 40000 K at least at some point for stars 
more massive than about 4OM0. Two example stellar birth lines from Bernasconi & 
Maeder (1996), crudely modified as described above, are shown by dashed-dotted lines 
in Figure |l^. It can be seen that they do indeed diverge from the ZAMS sufficiently to 
explain the appearance of a star with an age of 10^ yr at such low temperatures. 



7.3. A Young Cluster 

Fey et al. (1995) suggested that star ^^3 is a member of a young cluster. We are now 
in a good position to test this suggestion, as we know the extinction to the nebula and 
have good colors for the stars. Unfortunately, we cannot perform a statistical test for 
an overdensity of stars close to the UC HII region because completeness is so complex 
and poorly determined in the vicinity of the nebula. 

The stars with H-K > 1 in Table | and Figure |a are #3, #4, #13, #14, and 



#24. These stars are marked in Figure 11. Since they are projected close to the UC 



HII region and have extinctions that are similar the UC HII region, these are good 
candidate members of a cluster that is physically associated with the UC HII region. 
It is worth noting that with the exception of #3 and #14 these stars are not especially 
prominent in the images. By a visual inspection of the images alone, we might well 
have suggested that the bright stars to the NNE of star #3 were part of the cluster, 
but these appear to be foreground stars by virtue of their lower extinctions. 

The star with the most extreme infrared excess is star #4, located only 1.9 arcsec SSW 
of the ionizing star. This star is difficult to see without suppressing the nebula emission 
(see Figure ^) and has not been previously noted. Comparing it directly to the ionizing 
star gives an H— K excess of 0.75 with a 3a error of 0.19. The presense of this infrared 
excess star and the absence of similar stars away from the nebula is further evidence 
both that the cluster is real and that it is young. 

We estimated the spectral types of these stars from their H magnitudes as the K 
magnitudes of the infrared excess stars are presumably severely contaminated with 
emission from circumstellar dust. We scale the extinction to the nebula by a A^^ *^ law 
and find Ajj ~ 3.6. At a distance of 5 kpc the ionizing star #3 has Mh ~ —5.1 and 
the other stars #4, #13, #14, and #24 are much fainter with Mh ~ -3.6, -3.0, -2.7, 
and —1.8 with uncertainties of at least half a magnitude. Cluster stars other than 
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the ionizing star have luminosities appropriate for early to mid B stars (Schmidt-Kaler 

1982; Koornneef 1983). The situation in the chistcr appears somewhat similar to that 
in M17 (Hanson & Conti 1995) in which naked O stars appear with enshrouded B 
stars, although possibly here the division between naked and enshrouded stars is not 
so sharp. 

No stars are visible in the vicinity of the dense core located about 5 arcsec W of the 
ionizing star. This is consistent with the estimate of Ak ~ 400 through the core 
(Cesaroni et al. 1994). 

8. Summary and Discussion 

We have demonstrated that the stellar and nebular properties of UC HIT regions can 
be studied in the near infrared using moderate integration times on moderate aperture 
telescopes. These observations give direct, quantitative information on the color and 
luminosity of the ionizing star. Additionally, they give important information on the 
extinction to the nebula and the pattern of dust in and around the nebula. 

Our observations of G29.96— 0.02 confirm that the extended radio continuum emission 
seen by Fey et al. (1995) is indeed associated with the bright arc of emission seen by 
Wood & Churchwell (1989a). This extended emission does not appear to be naturally 
explained by the bow shock model. Fey et al. (1995) suggest a stationary model in 
which the arc and the extended emission are explained by the density gradient in the 
molecular cloud. The extinction to the UC HII region derived from radio continuum 
and Br7 imaging amounts to Ak ^ 2.25 or Ay ~ 25 and is higher in the region 
of the arc. The most likely explanation for the additional extinction is that we are 
seeing a dense molecular sheath around the arc or that there is significant internal 
extinction. The required magnitude of the internal extinction requires roughly twice as 
much column density in low density or neutral gas as there is in dense ionized gas in the 
vicinity of the arc. This might naturally occur in models in which UC HII regions ablate 
dense inclusions (Dyson 1994; Lizano & Canto 1995; Dyson, Williams, & Redman 1995; 
Redman, Williams, & Dyson 1996; Williams, Dyson, & Redman 1996; Lizano et al. 
1996) or in which the arc arises from a photoevaporation front being driven into an 
inhomogeneous medium. 

We confirm the identity of the ionizing star responsible for ionizing the UC HII region. 
It is located at the center of the arc of emission, as predicted by both the bow shock 
and stationary models for the morphology of the UC HII region. Direct light from the 
ionizing star accounts for only about 1/6 of the total light in each of the J, H, and 
K bands. The rest is presumably the result of thermal emission from dust, scattering 
from dust, and line and continuum emission from ionized gas. 
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The ionizing star does not possess a significant excess at K, indicating that it is no longer 
accreting. This also suggests that the photoevaporating disk model of Hollcnbach et 
al. (1994) may not be applicable to G29.96— 0.02 in its cmrcnt state, although detailed 
predictions of the spectrum of such a disk are required to confirm this. 

We confirm the suggestion of Fey et al. (1995) that the ionizing star is closely associated 
with a number of stars with luminosities appropriate for B stars. One of these stars 
has a strong infrared excess, confirming that the cluster is young. The lack of other 
luminous sources suggests that the ionizing star dominates the bolometric luminosity 
of the cluster and so is likely to have a mass of 60Mq or more. 

Our measurement of the intrinsic rriK of the ionizing star, along with existing limits on 
the bolometric luminosity and distance, allow us to place limits on where the ionizing 
star can appear in the HR diagram. We find that the ionizing star is too cool to be 
within 10^ yr of the ZAMS, in apparent contradiction with the age of about 10^ yr 
estimated for the the UC HII region. This suggests that the stellar birth line (Beech & 
Mitalas 1994; Bernasconi Ki Maeder 1996) for G29.96— 0.02 must be fairly cool: cooler 
than about 40000 K for some mass in excess of 4OM0. This places an important limit 
on how blue the stellar birth line can be; a red limit already exists, as the stellar 
birth line must extend to at least about 90Mq (Bernasconi & Maeder 1996). The 
implications of such a red birth line for the accretion history cannot be judged at 
the moment as theoretical models have not yet been constructed at the twice solar 
metallicity appropriate for G29.96— 0.02. We place a firm lower limit on the mass of 
the ionizing star at SOM© but favor a mass in excess of about 6OM0. 

Our discovery that the ionizing star in G29.96— 0.02 is significantly cooler than the 
ZAMS is a further strike against the method of determining the spectral type of the 
ionizing star in an UC HII region from an estimate of its ionizing continuum flux 
or spectral hardness. We can see this by considering the solar metallicity models of 
Schaercr & de Koter (1996). These models combine a stellar atmosphere code and 
a stellar evolution code. Their evolved 85Mq and 6OM0 models with logTeff ~ 4.6 
(models E3 and D3) have similar Lyman continuum fluxes to their ZAMS 4OM0 
and 25M0 models (models CI and Bl). Although these speciflc models have solar 
metallicity, we expect that similar discrepancies will arise regardless of metallicity. 
Attempts to use the bolometric luminosity to determine the mass of the ionizing stars 
in UC HII regions are in better shape, as massive stars evolve initially at almost constant 
Lboi, although source confusion remains a problem. 

This work has benefited from discussions with and comments by Paulo Bernasconi, 
Pepe Franco, Jay Gallagher, Margaret Hanson, Lynne Hillenbrand, Melvin Hoare, 
Todd Hunter, Tom Megeath, & Robin Williams and from a very thorough report 
from an anonymous referee. We thank Jon Holtzman for access to and instruction 
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in the use of his simultaneous PSF-fitting photometry routines. Our observations were 
made using the Perkins Reflector of the Ohio Wesleyan University and the Ohio State 
University at the Lowell Observatory We thank Andrew Afllerbach, Riccardo Cesaroni, 
Mark Claussen, Alan Fey, Ralph Gaume, and Stan Kurtz for allowing us to use their 
published and unpublished data. AMW and PH are grateful for partial support from 
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Figure Captions 



Fig. 1. — A 60 X 60 arcsec region around G29.96— 0.02. The axes are marked in arcsec 
north and west of the ionizing star. Darker shades indicate brighter emission, (a) J 
band, (b) H band, (c) K band, (d) Br7 with a white cross marking the location of the 
ionizing star. 

Fig. 2. — The spectral energy distribution for G29.96— 0.02. The points are taken from 
Table |2[ The dashed line is a u^'^ extrapolation from the 2 cm measurement and shows 
that optically thin free-free emission dominates to about 2 mm. The dotted line is a 
u'^ extrapolation from the 1300 /im measurement and shows that the sub-mm region is 
dominated by optically thin thermal dust emission. The solid line below 100 //m is a fit 
to v'^B^{T) which gives T = 26 K. The solid line above 100 /^m is a linear interpolation 
in (log F,,, log i^). 

Fig. 3. — A 30 X 30 arcsec region around G29.96— 0.02. In each panel, the white cross 
marks the location of the ionizing star, the axes are marked in arcsec north and west of 
the star, and the beam size is marked in the lower left. Darker shades indicate higher 
temperatures, (a) The brightness temperature Tb at 2 cm smoothed to the resolution 
of the 6 cm image. The contours are at 25, 50, 100, 200, 400, and 800 K. (b) The 
brightness temperature Tt, at 6 cm. The contours are at 250, 500, 1000, 2000, and 
4000 K. (c) The electron temperature Tg derived from the 2 cm and 6 cm images. The 
contours are at 4000, 5000, 6000, and 7000 K. (d) The brightness temperature Tb at 
2 cm smoothed to the resolution of the 21 cm image. The contours are at 25, 50, 100, 
200, 400, and 800 K. (e) The brightness temperature Tb at 21 cm. The contours are at 
500, 1000, 2000, and 4000 K. (c) The electron temperature Te derived from the 2 cm 
and 21 cm images. The contours are at 4000, 5000, 6000, and 7000 K. 

Fig. 4. — A 30 X 30 arcsec region around 029.96— 0.02. In each panel, the white cross 
marks the location of the ionizing star, the axes are marked in arcsec north and west of 
the star, and the beam size is marked in the lower left. Darker shades indicate higher 
values, (a) The observed flux in Br7. The contours are spaced by factors of two. (b) 
The observed flux in 2 cm radio continuum at full resolution. The contours are spaced 
by factors of two. (c) The observed flux in 2 cm radio continuum smoothed to match 
the resolution of the Br7 image. The contours are spaced by factors of two. (d) The 
apparent extinction at Br7. The contours are at 1.8, 2.0, 2.2, 2.4, and 2.6. 

Fig. 5. — A 30 X 30 arcsec region around G29.96— 0.02. In each panel, the white cross 
marks the location of the ionizing star, the axes are marked in arcsec north and west of 
the star, and the beam size is marked in the lower left. Darker shades indicate higher 
values, (a) The observed flux at 1.3 cm. The contours are spaced by factors of two. 
(b) The predicted flux at 1.3 cm from the 2 cm image. The contours are spaced by 



factors of two and have the same values as in (a), (f) The ratio of the observed flux to 
the predicted flux at 1.3 cm. The contours are at 0.8, 0.85, and 0.9 

Fig. 6. — A 30 X 30 arcsec region around G29.96— 0.02. In each panel, the white cross 
marks the location of the ionizing star, the axes are marked in arcsec north and west of 
the star, and the beam size is marked in the lower left. Darker shades indicate higher 
values, (a) The total line flux in C^^O. The contours are spaced by factors of two. 
Negative contours are dashed. The white circle marks the FWHM 4 arcsec ammonia 
clump seen by Cesaroni et al. (1994). The white crosses (x) mark the locations of the 
water masers seen by Hofner & Churchwell (1993). The white square (□) marks the 
locations of the formaldehyde maser seen by Pratap, Menten, & Snyder (1994). (b) 
The apparent extinction at Br7. The contours are at 1.8, 2.0, 2.2, 2.4, and 2.6. The 
C^^O contours from (a) are superposed. 

Fig. 7. — A 30 X 30 arcsec region around G29.96— 0.02. In each panel, the white cross 
marks the location of the ionizing star, the axes are marked in arcsec north and west of 
the star, and the beam size is marked in the lower left. Darker shades indicate higher 
values, (a) The model internal extinction at Br7. The contours are at 0.05, 0.10, 0.15, 
0.20, and 0.25. (b) The apparent external extinction after removing the model internal 
extinction. The contours are at 1.8, 2.0, 2.2, 2.4, and 2.6. (c) The apparent external 
extinction after removing twice times the model internal extinction. The contours are 
at 1.8, 2.0, 2.2, 2.4, and 2.6. 

Fig. 8. — A 60 X 60 arcsec region around G29.96— 0.02. Darker shades indicate brighter 
emission, (a) J band, (b) J band with nebular emission suppressed, (c) i^band. (d) H 
band with nebular emission suppressed, (e) Khand. (f) Khand with nebular emission 
suppressed. 

Fig. 9. — (a) Color-color diagram for stars within a 30 arcsec square centered on the 
ionizing star, (b) Color-color diagrams for stars outside that box. The error bars are 
la. Only stars with la errors of less that 0.15 in K and 0.3 in both colors are shown. 
The ionizing star is marked with a solid symbol. The dotted lines are the loci of main- 
sequence and giant stars. The dashed lines are reddening vectors of A~^-^, A~^-^, and 
^-2.0 reddening laws and correspond to Abtj = 2.20. The quadrilateral shows the 
expected colors of a hot star with Abtj = 2.20 ± 0.25 under these reddening laws. 

Fig. 10. — A theoretical HR diagram. The region allowed by la, 2a, and 3a 
uncertainties on the measured properties of the ionizing star are shown by thick solid 
lines. The thin dotted lines show the loci where my = 7.32 at different distances. Also 
shown are tracks and isochrones of the Z = 2Zq models of Meynet et al. (1994). The 
solid lines are tracks for 120, 85, 60, 40, 25, and 20Mq models, the dashed lines are the 
zero-age mam-sequence (ZAMS), and 1, 2, 3, 4 x 10^ yr isochrones. The dashed-dotted 



lines are approximate stellar birth lines. As described in the text, these birthlines are 
almost certainly somewhat too blue. 

Fig. 11. — (a) A 30 X 30 arcsec region around G29.96— 0.02 showing the K band 
emission. Darker shades indicate brighter emission. Candidate members of a cluster 
associated with the UC HII region are marked, (b) As (a) but with nebular emission 
suppressed. 
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Table Captions 

Table 1: Radio Data 

Table 2: Spectral Energy Didstribution 

Table 3: Stellar Photometry 



Table 1. Radio Data 



Authors 


Frequency 


Telescope 


Beam 


LAS 








arcsec 


arcsec 


Fey et al. (1995) 


14.9 GHz/2 cm 


VLA B+C+D 


0.6 X 0.5 


90 


Cesaroni et al. (1994) 


24.1 GHz/1.3 cm 


VLA C+D 


2.7 X 2.1 


60 


Afflerbach et al. (1994) 


4.87 GHz/6 cm 


VLA B 


1.8 X 1.5 


40 


Claussen & Hofner (1996) 


1.43 GHz/21 cm 


VLA A 


1.3 X 1.5 


40 


Hofner et al. (1996) 


110 GHz/2.7 mm 


OVRO 


3.8 X 5.4 


20 



Table 2. Spectral Energy Distribution 



Wavelength 




Reference 




Jy 




21 cm 


0.94 


Claussen & Hofner (1996) 


6 cm 


3.18 


Afflerbach et al. (1994) 


2 cm 


3.9 


Fey et al. (1995) 


1.3 cm 


3.53 


Cesaroni et al. (1994) 


7 mm 


2.85 


Wood et al. (1988) 


2.7 mm 


2.37 


Hofner et al. (1997) 


1.3 mm 


15.6^ 


Mooney et al. (1995) 


1.3 mm 


8.0 


Chmi et al. (1986) 


1.3 mm 


5.81*^ 


Mooncy ct al. (1995) 


790 ^m 


5.1 X 10^ 


Hunter (1997) 


350 ^m 


3.0 X 10^ 


Hunter (1997) 


100 iim 


1.15 X 10"^ 


IRAS PSC 


f\C\ i/rn 

\J\J Lt/L±l. 


7 50 X 10^ 


IRAS PSC 


25 /xm 


1.70 X 10^ 


IRAS PSC 


13.0 //m 


2.13 X 10^ 


Ball et al. (1996) 


12 /im 


2.17 X 10^ 


IRAS PSC 


12.6 fim 


1.84 X 10^ 


Ball et al. (1996) 


12.0 fim 


1.57 X 10^ 


Ball et al. (1996) 


9.7 /xm 


4.84 X 10^ 


Ball et al. (1996) 


8.5 /xm 


4.30 X 10^ 


Ball et al. (1996) 


5.0 /tm 


5.0 


Herter et al. (1981) 


3.8 /xm 


1.9 


Hertcr at al. (1981) 


2.2 /tm 


4.40 X 10-1 


This work 


1.65 /tm 


9.48 X 10-2 


This work 


1.25 /tm 


1.66 X 10-2 


This work 



^Integrated flux 

^12 arcsec FWHM point source 



Table 3. Stellar Photometry 



Star 


Ax 
arcsec 


Ay 
arcsec 


J-H 


H-K 


K 


1 




9'^ Q 


_i_n 87 -1- n ni^ 


-rU.4tl It u.uo 


y.ui It u.u^ 


9 


-U07 A 


OA fi 


If) Qfi -I- n nfi 


T^u.4/ It u.uo 


lU.lZ It U.U4i 


o 
o 


-ru.u 


-pU.U 


1 9 /in -1- n ni^ 
-pz.'iU It u.uo 


-rJ-.oy It u.uo 


iU.oO It U.U'i 


4 


+1.6 


— 1.0 


1 O O O 1 Mil 

+2.28 ± 0.11 


+2.34 ± 0.07 


11.14 ± 0.05 


5 


+13.8 


-3.4 


+1,22 + 0.05 


+0.57 + 0.05 


11.41 + 0.03 


6 


-18.8 


-19.5 


+1.68 + 0.05 


+0.89 + 0.04 


11.55 + 0.03 


7 


+3.9 


+4.8 


+1.84 + 0.05 


+0.83 + 0.05 


12.34 ± 0.04 


8 


+5.1 


+9.5 


+0.76 ± 0.05 


+0.44 ± 0.05 


12.35 ± 0.04 


9 


+3.2 


-22.3 


+0.63 ± 0.04 


+0.37 ± 0.05 


12.46 ± 0.03 


10 


+2.5 


+14.3 


+1.89 + 0.06 


+0.86 + 0.06 


12.92 + 0.05 


11 


-25.6 


+8.7 


+1.13 + 0.05 


+0.74 + 0.05 


12.93 + 0.03 


12 


-1.4 


-24.0 


+2.02 ± 0.08 


+0.95 + 0.06 


13.10 + 0.04 


13 


-5.5 


+13.8 


+2.79 + 0.26 


+1.86 + 0.09 


13.20 ± 0.05 


14 


-4.4 


+3.7 


+1.80 ± 0.09 


+1.15 + 0.08 


13.23 ± 0.06 


15 


-3.8 


-4.8 


+0.98 ± 0.06 


+0.79 ± 0.08 


13.41 ± 0.06 


16 


+23.4 


+6.7 


+2.04 + 0.07 


+0.97 + 0.06 


13.42 ± 0.05 


17 


-22.9 


-6.3 


+1.69 + 0.06 


+0.83 + 0.06 


13.47 + 0.05 


18 


-2.4 


+1.1 


+1.34 + 0.15 


+1.66 + 0.14 


13.56 + 0.11 


19 


+16.2 


+23.9 


+0.72 ± 0.06 


+0.44 ± 0.07 


13.61 ± 0.05 


20 


-2.3 


-2.9 


+1.92 + 0.22 


+1.80 + 0.17 


13.68 + 0.13 



Table 3 — Continued 



Star 


Ax 
arcsec 


Ay 
arcsec 


J-H 


H-K 


K 


21 


+17.6 


-1.1 


+2.28 ± 0.38 


+2.18 + 0.21 


13.68 ± 0.09 


22 


+18.8 


+7.2 


+2.84 + 0.41 


+2.07 + 0.11 


13.75 + 0.08 


23 


-0.6 


+6.1 


+0.70 ± 0.05 


+0.23 + 0.09 


13.79 + 0.08 


24 


-10.0 


-3.2 


+1.97 + 0.17 


+1.39 + 0.13 


13.93 + 0.10 


25 


+22.7 


-3.3 


+0.52 ± 0.06 


+0.08 ± 0.12 


14.01 + 0.11 


26 


-7.8 


+0.1 


+1.04 + 0.07 


+0.62 + 0.11 


14.03 ± 0.10 


27 


-29.3 


+13.4 


+1.38 + 0.08 


+0.83 ± 0.07 


14.06 ± 0.05 


28 


-17.6 


-1.9 


+1.14 + 0.40 


+2.90 + 0.27 


14.25 ± 0.08 


29 


+ 15.3 


+ 19.2 


+1.58 + 0.33 


+1.84 + 0.20 


14.36 + 0.08 


30 


+29.5 


+28.5 


+0.95 ± 0.09 


+0.48 ± 0.09 


14.41 ± 0.07 


31 


+8.4 


+5.1 


+0.79 ± 0.06 


-0.03 ± 0.17 


14.48 ± 0.16 


32 


-9.4 


+6.2 


+2.25 ± 0.22 


+0.95 ± 0.16 


14.49 ± 0.14 


33 


+3.5 


+29.0 


+1.01 + 0.11 


+0.82 + 0.12 


14.49 ± 0.10 


34 


-13.2 


+17.5 


+1.11 + 0.09 


+0.42 + 0.12 


14.53 + 0.10 


35 


+20.0 


+6.0 


+1.07 + 0.10 


+1.13 + 0.14 


14.56 + 0.12 


36 


-20.3 


-13.8 


+1.89 + 0.26 


+1.23 + 0.15 


14.83 + 0.10 


37 


+29.1 


+1.3 


+1.36 + 0.11 


+0.72 ± 0.11 


14.83 ± 0.08 


38 


-26.4 


+23.2 


+2.36 ± 0.53 


+1.43 + 0.18 


14.87 + 0.11 


39 


-27.3 


+10.1 


+1.58 + 0.19 


+0.85 + 0.16 


14.91 + 0.12 


40 


+0.7 


-22.0 


+1.48 + 0.26 


+1.06 + 0.18 


14.94 + 0.12 



Table 3 — Continued 



Star 


Ax 
arcsec 


Ay 
arcsec 


J-H 


H-K 


K 


41 


-12.4 


-17.3 


+0.74 ± 0.10 


+0.32 ± 0.13 


14.98 ± 0.10 


42 


-25.9 


-9.3 


+1.23 + 0.15 


+0.73 + 0.16 


15.01 + 0.13 


43 


+19.1 


+14.4 


+0.64 + 0.10 


+0.26 + 0.16 


15.06 + 0.13 


44 


-0.6 


+14.6 


+0.90 + 0.10 


+0.03 + 0.29 


15.24 + 0.28 


45 


-20.5 


-19.7 


+0.70 + 0.12 


+0.18 + 0.24 


15.34 + 0.22 


46 


-7.6 


-27.4 


+0.62 + 0.12 


+0.28 ± 0.19 


15.42 + 0.17 


47 


-17.6 


+1.1 


+1.15 + 0.11 


+0.34 ± 0.26 


15.52 ± 0.25 


48 


+14.2 


-14.5 


+1.20 + 0.12 


+0.35 + 0.20 


15.53 + 0.18 


49 


+3.9 


-6.2 


+2.10 + 0.34 


+0.48 ± 0.43 


15.60 + 0.41 


50 


-21.5 


-17.9 


+1.61 + 0.45 


+1.22 + 0.28 


15.61 + 0.19 


51 


-25.6 


+17.8 


+0.70 + 0.11 


+0.14 + 0.20 


15.65 + 0.17 



